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Storing quantum information for long times without disruptions is a major requirement for most 
quantum information technologies. A very appealing approach is to use self- correcting Hamiltonians, 
i.e. tailoring local interactions among the qubits such that when the system is weakly coupled to a 
cold bath the thermalization process takes a long time. Here we propose an alternative but more 
powerful approach in which the coupling to a bath is engineered, so that dissipation protects the 
encoded qubit against more general kinds of errors. We show that the method can be implemented 
locally in four dimensional lattice geometries by means of a toric code, and propose a simple 2D 
set-up for proof of principle experiments. 



I. INTRODUCTION 

There are two existing approaches to providing coher- 
ent quantum storage on many-body systems. The first 
one corresponds to fault tolerant quantum circuits [TJ [2] . 
If one can perform quantum gates and provide fresh 
initialized qubits with a sufficiently high accuracy and 
frequency, then quantum computing and in particular, 
quantum memory is possible for a time exponential in 
the dedicated resources. 

More recently, Kitaev [3J H] proposed that it might 
be possible to protect quantum information passively 
by engineering of suitable Hamiltonian systems, in anal- 
ogy to magnetic domains for classical memories. While 
an energetically degenerate code subspace insensitive to 
Hamiltonian perturbations is a necessary condition, it 
has become clear that there are additional requirements 
for this approach to quantum memories to work. Possi- 
bly the most important requirement is to cope with the 
undesired coupling between the storage system and its 
environment. In this direction, the approach that has 
benefited from the most theoretical progress goes by the 
moniker of self-correcting Hamiltonians, [SHE]- 

For self-correcting Hamiltonians, a weak local coupling 
to a thermal bath is assumed. Making a Born-Markov 
approximation, the evolution of the system can be de- 
scribed by a thermalizing master equation. While for 
general local couplings, any initial state will decay to the 
unique Gibbs state, it is still possible for the decay rate of 
specific observables to become smaller as the number N 
of subsystems increases. This leads to the possibility of 
storing quantum information by encoding it on a pair of 
slowly decaying anticommuting many-body observables. 
A Hamiltonian will thus be called self-correcting pro- 
vided that below a certain finite bath temperature the 
dissipative dynamics leads to information lifetimes grow- 
ing with the system size (typically following an expo- 
nential increase). Alicky et al. [6] rigorously proved an 
exponentially long relaxation time for protected observ- 
ables in the 4D toric code. Chesi et al. [H] generalized 
this result deriving a criteria for quantum memory based 
on self-correcting Hamiltonians and lower bounds on the 
storage times. However, it is in general not known how 
non thermal noise or even thermalization under a per- 



turbed Hamiltonian [TO] affects this lifetime. In particu- 
lar, this may be the case whenever the qubits are weakly 
coupled to an additional bath which induces a small rate 
of depolarization[TTj . 

Building on previous results, we propose and analyze 
an alternative way of protecting quantum states. The 
method is similar to that of protecting Hamiltonians, but 
now the main idea is to tailor the coupling of the qubits 
to a bath, so that the engineered dissipation extends the 
life-time of the encoded qubit. Apart from being passive 
(i.e. not requiring the burden of interrogating the quan- 
tum memory at intermediate times) , the main advantage 
of this scheme is that it can potentially correct for other 
kinds of errors beyond those generated by thermalization, 
including depolarizing noise. In particular, we propose a 
specific method in 4 spatial dimensions inspired by toric 
codes and obtain evidence of its performance with the 
help of numerical simulations. We also investigate a sim- 
plified 2-dimcnsional model protecting only from phase 
errors which could be a good candidate for proof of prin- 
ciple experiments. 

Many-body classical memories based on dissipation 
(often under the name of asynchronous celular automata) 
have naturally appeared in the context of classical fault 
tolerant computation. For example, using a simple lo- 
cal update rule on a 2D lattice, Toom [T3] showed 
that classical information can be protected against weak 
local noise. A more elaborate update rule by Gacs [T3] 
provide protection even on a ID lattice. These results 
already suggest that dissipation may offer a powerful al- 
ternative to the existing methods for constructing many- 
body quantum memories, as investigated in the present 
work. In fact, several authors have already proposed the 
use of continuous quantum error correcting codes [ToTl"T§] . 
However previous works concentrate on a single level of 
error correction and do not address the large N many- 
body scenario. A notable exception is the work of Dennis 
et al. [3J introducing a heat bath algorithm (thermal dis- 
sipation for the 4D toric code) in order to simplify the 
efficacy analysis of a local many body quantum error cor- 
rection algorithm. At the crux of this approach is that 
thermal dissipation can be interpreted not only as in- 
troducing decoherence (errors), but also as performing a 
form of error correction, with the balance between the 
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two effects roughly given by the bath temperature. In- 
deed, this heat bath algorithm can already be seen as a 
dissipative quantum memory lending itself to more nat- 
ural engineering. In fact, engineered dissipation is more 
general in that it need not satisfy detailed balance con- 
ditions and thus its power extends that of cooling a self- 
correcting Hamiltonian. In other words, the steady state 
need not be an equilibrium state and its dynamics may 
show a net flow (imagine a funnel receiving water from 
a hose). As the classical results show, this more gen- 
eral kind of dissipation may be crucial in order to correct 
general kind of errors. 

Our proposal can be viewed as another example where 
engineered dissipation may become a useful and alterna- 
tive tool in the context of quantum information process- 
ing, beyond quantum computation , state engineering 
[2Ql[21], or entanglement creation [22 . In all those cases, 
it is desirable to be able to couple small subsets of qubits 
to Markovian environments so that their evolution equa- 
tion follows a prescribed master equation. As exposed 
in |20j . dissipative gadgets provide a direct way of im- 
plementing this is in terms of damped qubits; that is, a 
set of qubits which themselves follow a damping master 
equation due to their coupling to an environment. Those 
qubits can be directly coupled to the physical qubits of 
the quantum memory or computer to provide the desired 
dissipation, and thus appear as an important resource in 
dissipative quantum information processing. 

This paper is organized as follows. In Section II we 
briefly present the general idea of engineered dissipative 
quantum memories. In Section III we display two differ- 
ent but rather obvious approaches to dissipative quantum 
memories and discuss why they are not entirely satisfac- 
tory. In Section IV we present a specific method in 4 
spatial dimensions as well as the results of numerical sim- 
ulations which validate the performance of the scheme. 
Section V contains a simplified version in 2 spatial dimen- 
sions which corrects against phase errors and that could 
be tested experimentally in the near future. In Section 
VI we show how one can use dissipative qubits to engi- 
neer the dissipation and analyze under which condition 
one can use them in this context. All previous section 
contain the main statements of our work. The detailed 
proofs of our results and more thorough explanations are 
given in the appendices. 



II. STATEMENT OF THE PROBLEM 

We consider a logical qubit encoded in N physical 
qubits, which are appropriately coupled to an environ- 
ment providing dissipation. We describe the action of the 
engineered environment, as well as of the other sources 
of decoherence through a master equation 



P = £diss(p) + £noise(p) 



(1) 



Here, p is the density operator for the qubits, £diss the Li- 
ouvillian describing the engineered dissipation, and £ n oise 
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FIG. 1. (color online) We assume that a piece of quantum 
information is encoded into a many body system. The engi- 
neered dissipation, is then responsible for making the degrees 
of freedom which carry the encoded quantum information re- 
silient against the uncontrolled noise processes taking place. 
Finally, the decoding process extracts the quantum informa- 
tion from the collective degrees of freedom. 



will denote a noise term contribution to the master equa- 
tion. This could be local depolarizing noise for instance 



£noise(p) = ^e^dep(p) — E 



N 1 

T 

n=l 



itr„(p)-p, (2) 



or any other weak local noise term. Our goal is to show 
that for appropriate choices of £diss we can extract the 
encoded qubit reliably after a time which substantially 
increases with N. 

In general, any trace preserving dissipative master 
equation as £diss may be write in Lindblad form |23j 

p = C(p) = -i{H,p] + £ 2L kPa L\ - {L\L kP } + , (3) 



consisting of a Hamiltonian term describing the unitary 
evolution, and a dissipative part which may be written in 
terms of Lindblad or jump operators L k . Furthermore, 
the models of engineered dissipation we propose can be 
seen to adhere to a more benign form 



(4) 



where Ti are positive trace preserving channels. For these 
particular cases, the time dependent density matrix may 
be given an explicit stochastic expansion in the form of 



Pit) 



E 



T n p(0) 
n! 



(5) 



where r = J2i r * an d T(p) = J2i ^iTi(p)- This stochastic 
expansion will be useful for both proofs and Monte Carlo 
simulations. 



III. STRAIGHTFORWARD QECC ENCODING 

Here we introduce and analyze two straightforward 
methods of encoding a QECC in the dissipation. The 
first one consist of coupling all the qubits with a reservoir 
in such a way that each application of a jump operator 
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a whole error correction procedure takes place. In the 
second, we encode the QECC in several Lindblad terms, 
so that each jump correspond to an execution of a part 
of the QEC. The main purpose of this section is to show 
that those simple approaches do not work as one could 
imagine, and thus it illustrates why the design of engi- 
neered quantum memories is not a trivial task. Both 
approaches require multibody coupling to a single en- 
vironment, where the number of system qubits coupled 
to the same damped qubit grows with N, the size of the 
memory. In principle perturbation theory gadgets allow 
the engineering of such terms, provided their respective 
intensity decay exponentially with the number of subsys- 
tems involved. Not withstanding, a strength increasing 
with N would be required to make the first approach 
work, while in the second approach only a polynomial 
decrease with the number of subsystems involved would 
preserve functionality. In the next section we will present 
a scheme which circumvents these problems, although 
still with the caveat that it requires non-local couplings 
(as it works in 4 spatial dimensions). 



A. Single Jump Operator 

One major obstacle to traducing the usual error cor- 
rection strategies to a dissipative scenario is due to the 
random times at which dissipative terms enact the re- 
covery operations. We illustrate this problem in the case 
of a straightforward approach to dissipative protection. 
One can always implement in the dissipative Liouvillian a 
standard quantum error correction procedure which pre- 
serves the logical qubit: Cd lss {p) — T[R.(p) — p], where 
1Z is a full recovery operation and T adjusts the rate at 
which the recovery operation is applied (imagine full cor- 
rection of an N qubit QECC) . Apart from the unrealistic 
nature of highly many-body dissipation terms required 
in this construction, it is easy to see that it does not 
serve our purposes. The reason can be seen by unravel- 
ing the quantum jump operators [23], there is a finite N 
independent probability for more than i- time to elapse 
until the next recovery operation. Such long times allow 
too many errors to accumulate for any QECC to recover 
with high fidelity. 

The alternative is to have dissipation implement many 
independent processes instead of a single monolithic er- 
ror correction procedure. Ideally, having independent 
processes take care of removing independent error sets 
can make the accumulation of a critical fraction of errors 
exponentially unlikely. The difficulty of having indepen- 
dent dissipation processes is that contrary to the circuit 
model the order of their application is not enforced in 
any way. Thus, directly encoding each gate of a QECC 
recovery circuit into a dissipation term generally leads to 
a meaningless evolution. However, we will show that in 
specific cases where dissipation terms commute or show 
some order property lending itself to rigorous analysis, 
the asynchronous nature is not an obstacle. 



B. Concatenated QECC Dissipation 

It is indeed possible to design a many-body dissipative 
quantum memory. The strategy is to take the dissipa- 
tion term as a sum of recovery operations occurring on 
different groups of qubits. Those operations correspond 
to recovery of the different logical qubits at each level of 
a simple concatenated QECC [35] ■ Intuitively, one may 
argue that the difficulty of implementing a given dissipa- 
tion term increases with the number of qubits involved. 
We attempt to compensate for this difficulty by impos- 
ing that the operator norm required for such Lindblad 
terms decays with a power law respect to the number of 
physical qubits involved. More specifically, we take 



£ diss ( P ) = rJ2s M - l [n ltn {p)-p\. 



(6) 



l:n 



Here, I — 0, 1, . . . , M — 1 denotes the level of concate- 
nation, and n further specifies on which set of qubits 
the recovery operations lZi. n are applied. In appendix 
|X) we show that if the local noise rate r e is sufficiently 
small then initially encoded information is lost at a rate 
which is exponentially small with respect to the number 
of qubits used (i.e. double exponentially small with the 
level of concatenation M). The weakness condition on 
the noise can be made precise by 



r £ <r; 



k 2 ' 



(7) 



where k is the number of physical qubits in the code to be 
concatenated. Assuming the perfect 5 qubit QECC and 
taking the strength of many body terms inversely pro- 
portional to the number of bodies (5 — 1/5), a threshold 
of r* = 1.6 x 10~ 3 r is obtained for the noise rate. When 
the error rate is below the error threshold, the relaxation 
rate for the encoded information has an exponentially 
small upper bound given by 



7-- 1 < T c 5 



(8) 



The above scheme is mainly of formal interest, since 
the non local recovery operations encoded in the dissi- 
pative master equation require many qubits at different 
locations to interact with the same environment. While 
the necessary scaling of such terms needs to be polyno- 
mial for our proof to go through, the derivation of such 
terms based on effective many-body Hamiltonians and 
the dissipative gadgets we propose is expected to decay 
exponentially with the number of bodies involved. Even 
more realistically, one would expect many-body dissipa- 
tion terms to cope with many body error terms arising 
from imperfect implementation. In practice, it would be 
desirable to find a set up where the dissipation terms are 
spatially localized by considering the qubits arranged in 
a lattice. 
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IV. LOCAL DISSIPATIVE PROTECTION IN 4D 

In classical systems Toom's rule |12j has been proven 
to be a simple translationally invariant update rule in 
a 2D Periodic Boundary Condition (PBC) lattice which 
is capable of preserving classical information, provided 
that the noise contribution to the dynamic is sufficiently 
weak. While we have not been able to extend this rule for 
quantum protection in 2D, we will consider a quantum 
analog of Toom's rule for 4D. The underlying QECC used 
is the 4D toric code, a stabilizer quantum error correct- 
ing code with 6 body stabilizer generators which can be 
made spatially local in a 4D PBC lattice. Dennis et al.[5] 
proposed it as a local QECC, and the corresponding sta- 
bilizer Hamiltonian was recently rigorously proven to be 
self- correcting by Alicki et al. [5] . We derive a local mas- 
ter equation for protecting information encoded into the 
4D toric code based on a Toom like rule introduced by 
Ahn [26] and study its efficiency for protecting encoded 
observables. We then consider the protection process and 
numerically study the lifetime of information when depo- 
larization errors are introduced extensively at a small yet 
constant rate. 

A fully rigorous description of the QECC and the local 



update rule used is provided in the appendix IX For the 
moment it is sufficient to specify that the master equation 
has the form of cq. where the specific £diss used 
associated to the 4D toric code will be called £4D-TCTooro 
and £ n oisc is weak extensive depolarizing noise as in eq. 
([2]). The numerical results (Fig. [2| strongly support 
the existence of a critical error rate T* « 0.004 x T ( 
where T is the correction rate to be specified ) below 
which, the lifetime of the encoded information increases 
exponentially with the lattice size. 

Although the results above have no obvious practical 
implication, they suggest that local models may exist in 
spatial dimensions lower than 4 (for the search of quan- 
tum memories based on protecting Hamiltonians in lower 
dimension see [FJ [7J [8]). The hope, is that even if self- 
correcting quantum memories fail to exist in lower dimen- 
sions, the use of engineered dissipation may still provide 
a solution. 



Numerical simulations 
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FIG. 2. (color online) The mean time to error for a logical 
observable is plotted in log scale units of p. Error rates T e 
are provided in units of V. The plots further suggests the ex- 
istence of a critical value for error rates F* w 0.004. (a) Each 
curve corresponds to a fixed odd value of the lattice size N. 
The independent axis F e /r is also in log scale suggesting that 
for each fixed N the information lifetime show an asymptotic 
(small T e ) power law dependence with 1/F e with the expo- 
nent increasing for larger N. (b) Each curve corresponds to a 
fixed value of the error rate IV For low error rates r e < F*, 
lifetime is seen to improve exponentially with N. 



The key feature that allows us to perform efficient sim- 
ulations of the relaxation times for logical observables, is 
that the terms in £4D-TCToom may be naturally split in 
two subsets, where terms in one subset commute with 
terms in the other. Thus, efficient classical Monte Carlo 
simulations provide unbiased estimates for expectation 
values and correlations for half of the stabilizers and 
half of the logical observables. Throughout each simu- 
lation the relevant error corrected logical observable was 
measured on a copy of the system state after every unit 
of simulated time. Simulation were interrupted when a 
measurement outcome differing from the initial value was 



obtained. For each parameter, lattice size N and the de- 
polarization rate T e , a total of 1000 such runs were per- 
formed and the simulated times were averaged to obtain 
the relaxation time presented. These simulations where 
performed on 62 AMD Opteron processors taking a total 
of five days to obtain the data presented (Fig. [2]). 

V. ACCESSIBLE TOY MODEL 

As a a proof of principle, we now present an engineered 
dissipation toy model providing protection for quantum 
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information. One can implement the underlying ideas 
of dissipative quantum memories with 2D lattices at the 
expense of being able to correct only for dephasing noise 
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(9) 



can increase the relaxation time of X EC , a many-body 
encoded observable (red) . This is in contrast to the com- 
plementary observable which does not benefit from dis- 
sipative protection. On the contrary, given any depolar- 
ization rate, the relaxation time of Z L decreases with the 
inverse of the number of physical qubits involved (blue) . 



Given a noise model including only one type of error 
(such as <j z phase errors) we will be able to cast a classi- 
cal memory prescription into a quantum scenario. A first 
step, is to define two logical observables 



zEC 



X 



EC 



(10) 



where 9 is the Heaviside step function. The first observ- 
able Z commutes with the noise £ p hase and is thus 
completely immune to it. The noise can only change 
the value of X EC , for the part of p which is in the ±1 
eigenspace of ^ s erf (i.e. states for which the absolute 
magnetization in the X direction is minimal). Dissipa- 
tion will protect the X EC observable by keeping most of 
p in a high X magnetization subspace. The master equa- 
tion p = £nn(p) for nearest neighbor majority voting is 
written as a Liouvillian in Lindblad form |23j as 

£nn(p) = T ^ L s ^xph\_ r t - tr>t L St r !t ,p} + , 

<s,r,t> 

(11) 

where the index s runs over all sites, r ^ t are nearest 
neighbors of s and the Lindblad operators are given by 



L, 



1 



i off 1 



(12) 



This is, the first factor performs a phase flip when the 
second and third factors (projectors) are non zero (i.e. 
when site s points differently than its two neighbors r 
and t). The Lindblad operators are designed such that 
they also commute with Z EC and can only change X EC 
in the portion of p with minimal X magnetization. 

The stability of the X EC observable in such an evo- 
lution can be mapped to magnetization metastability in 
classical studies c 27_, 28] . Restricting r and t to be north 
and east neighbors in an N x N PBC lattice, one recov- 
ers Toom's rule [121 113] which is proven to provide an 
exponential survival time, even in the presence of biased 
errors. However, the PBC requirement is experimentally 
unrealistic. 

We numerically consider an experimentally accessible 
setup which does not require periodic boundary condi- 
tions. Physical qubits will be located on an N x N 2D 
square lattice sites. The sites r and t are taken among 
all possible nearest neighbors of s. The number of valid 
neighbor combinations are ( 2 ) = 6 for inner sites s, 
( 2 ) = 3 for lattice border sites s and only one combi- 
nation for corner sites. In the following plot (Fig. |3|, 
we show how having a protective dissipation term Cnn 
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FIG. 3. (color online) Relaxation time for Z EC (red curves) 
and X L (blue curves) in units of F _1 . Each red curve presents 
the relaxation time tz (numerically obtained) corresponding 
to one value of the relative dephasing rate Y/Y phase given by 
the intercept at N — 1. Blue curve have the functional form 
tx = Y~* p * N~ 2 and each corresponds to one value of Y/Yd ep 
also given by the intercept at TV = 1. The lifetime r of the 
encoded logical qubit can be seen to be estimated by r ~ 
min{jx , tz} ■ Given Y/Ydep and Y/Y p hase, one may intersect 
the corresponding curves to obtain the value of N leading to 
the optimal qubit lifetime r. For example, if Ydep = 5 x 10 F 
and F phase = O.ir the optimal lattice size of 4 x 4 allows a 
x 100 increase in the quantum information relaxation time r. 
A more extreme case may be seen when Y p hase = O.Oir and 
Ydep < 5 x 10~ 5 r where a factor x50 is gained by simply 
using a 2 x 2 lattice. 



VI. DISSIPATIVE GADGETS 

As we have shown, the possibility of controlled quan- 
tum dissipation opens a host of new possibilities for QIP 
[S0Tl22j . However, while some naturally occurring forms 
of dissipation may be readily exploited, it is crucial to 
have a systematic way of engineering arbitrary dissipative 
dynamics. A way of achieving complete control over the 
dissipation is to be capable of engineering independent 
Lindblad jump operators while keeping their interference 
with each other weak. For this we must assume avail- 
ability of many body Hamiltonians, achievable through 
perturbation theory gadgets [29l |30] and of some natu- 
rally occurring dissipation, namely in the form of damped 
qubits. We apply the approximation of independent rates 
of variation [31] pg. 356 on the damped qubits which 
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requires the bath correlation time for the damping pro- 
cess to be much shorter than the inverse of any coupling 
constant in the system. Coupling to these damped qubits 
can thus be seen as a resource in the design of quan- 
tum dynamics, analogous to freshly initializing qubits in 
quantum circuits. 

Coupling the system to a damped qubit ancillary degree 
of freedom was proposed as a possible path to engineer 
arbitrary effective dissipative dynamics [20] ■ More specif- 
ically, the Hamiltonian coupling H = w{L®a + ®a~) 
to an ancilla with damping rate 7 leads to an effec- 
tive dissipative dynamics of the system corresponding to 
the Lindblad operator cuy/2/jL. Here a~ = |0) (1| and 
0+ = |1) <0|. 

In order to use these dissipation gadgets as basic build- 
ing blocks in more complex scenarios, it is essential to 
make explicit possible limitations and restrictions of the 
implemented dissipation. In appendix |VIII| we provide 
a detailed derivation of the effective system dynamics 
which makes three main contributions to our understand- 
ing of dissipative gadgets. Firstly, while the usual ap- 
proach of adiabatic elimination obtains an effective dy- 
namics in terms of a coarse grained time, our derivation 
shows that excluding a short initial transient period, this 
temporal coarse graining is not necessary. Secondly, we 
provide explicit bounds on the deviation from the desired 
state and instantaneous dynamics which are accompa- 
nied by a smallness prefactor (w/7) 2 . Finally, we in- 
clude an independent internal dynamic for the system, 
and show, that the resulting effective dissipation carries 
through essentially unaffected provided the strength of 
the internal dynamics is sufficiently weak. While the last 
point already suggests that the extensive application of 
local dissipation gadgets should be well behaved, a fully 
rigorous analysis is beyond the scope of this article. 

VII. CONCLUSIONS AND PERSPECTIVES 

We have introduced engineered dissipation as a tool 
to protect against general quantum noise and proposed 
examples providing protection from local noise. In the 
case of concatenated code dissipation, we prove that in- 
formation can be made resilient against any strictly lo- 
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VIII. ADIABATIC ELIMINATION OF ANCILLA 

In this section, we prove that Master equations with 
arbitrary Lindblad operators may be engineered to high 
accuracy by coupling the system to ancillary resource 
qubits which are themselves being cooled. The basic 
idea is to extend the system with an additional binary 
degree (spin 1/2) of freedom per Lindblad operator L to 
be implemented. These degrees of freedom are further 
assumed to be strongly dissipatively driven with a rate 
7 into a |0) (0| ground state. We will show that a target 
dissipative evolution composed of a single Lindblad jump 
operator 

Aargct(p) = LpL f - \ {tfL, P } + , (13) 

may be implemented within a small error margin. The 
technique used for the proof follows the adiabatic elimi- 
nation of the excited ancilla subspace in spirit, but takes 
into account corrections in order to provide rigorous 
bounds on the deviations from the intended evolution. 

Our derivation starts by assuming that the full dynam- 
ics of the system can be written as 

p = -i[H,p]+2 7 a-pa + - 7 {a + cr-,p} + + 1 C sys (p) (14) 

where a + = |1) (Oj^ and <j + = |1) (0| A are raising and 
lowering operators on the ancilla qubit and the Hamilto- 
nian H couples the system to the ancilla 

H = uj(L®(t + + ®cr-) (15) 

and £ sys is an additional evolution term with no effect 
on the ancillas. Here, the assumption that is implicitly 
being made, is that we may independently sum the in- 
teraction Hamiltonian H to the dissipative dynamics on 
both the system and the ancilla. In the case of the ancilla 
decay this is the approximation of independent rates of 
variation |31j pg. 356, which assumes correlation times 
for the reservoir responsible for spontaneous decay to be 
much shorter than any other relevant time in the sys- 
tem. An important example where this approximation 
holds to a great degree of accuracy is for two level atoms 
at optical frequencies, where the autocorrelation time of 
the coupled vacuum fluctuations can be as much as ten 
orders of magnitude shorter than the inverse of any of 
the other coupling constants. Since our derivation for 
the weak system Liouvillian does not require temporal 
coarse graining, the successively incorporation of Hamil- 
tonian interactions rigorously leads to the additive ap- 
pearance of the desired Liouville terms up to leading or- 
der. Assuming e — lo/j <C 1 we can rescale to a unitless 



time by incorporating a factor 7 leading to the following 
differential equations for the reduced density matrices. 

poo := (0| p |0) = 2/9ii - ie^Vio + itpoiL + £ sys (poo) 

(16) 

P01 := (0| p |1) = -poi + iepoo-k 1 ' - ietf p n + £ sys (poi) 

(17) 

pn := (1| p |1) = -2pn - ieLp i + iep 10 L^ + C sys (pn) 

(18) 



From here, we may obtain the integral forms 



such manipulation, one reaches the expression 
^2 r 



L Poo (t)L^ - e- T (2 - e- T )Lp 00 (0)tf (24) 
e~*'(2 - e-*')L/5 o(r - t')LUt' 



'{2-2e- t ')LL^p 11 {T~t')dt' 







h.c. 



le I -2t' 

2' e 



LC sys [p i(T - t')]dt' + h.c. 







-it' 



C ays [ Pll {T-t')}dt' 



Poi(t) =e T p O i(0) + J e 1 £ sys [p i(r - t')]dt' (19) Using eqs . ^ an d ([24^ one may prove tne following 



+ie / e- 1 ' [p m {r - t')ti - Pii {t - t')}dt' 



pn(r) =e- 2 >n(0)+ / e- 2r 'C sys [ Pll (r - t')]dt' (20) 
Jo 

-ie f e- 2t ' [L Pqi (t - t') - piq(t - t'JLtjdt' 



Assuming the initial conditions poi(O) = /On(0) = 0, 
that ||L|| < 1 and ||>C sys || < Ee 2 , and using that 
llpooll + HpuII < 1 we may bound 



IPoiWII < e and ||pii(r)|| < l 2 



(21) 



with e = 1 _% e 2 ■ It is now straightforward to bound 
||poo( r )ll < {A+E)e 2 . We may now concentrate on tighter- 
bounds composed of higher order terms in e but also, of 
exponentially decaying terms. A first step to do this is 
to perform integraton by parts; on eq. ( 19 1 one obtains 



poi (r) =zepoo (r)i f - iee~ r p 00 {0)L^ (22) 
*' [p m {T-t')L^ +L^p 11 {T-t , )]dt' 



-it I e 
'0 



+ e £ sys [p i(r - t )]dt . 
Jo 

In the case of p\\ we straightforwardly obtain 



Pn(r) 



: Lpoi(r) 



(23) 



2 Jo 



,-2t' 



Lpoi(r - t')dt' + h.c. 



-21' 



C sys [ Pll (T - t')}dt' 



higher order bounds 

||poi - iep m tf\\ < (2E + 5)? 3 + e e " T (25) 
\\ Pll - e 2 Lp m tf\\ < (3E + 7)? 4 + 2e 2 e~\ (26) 

Inserting these bounds into the definition of poo we may 
bound deviation from the target evolution by 

||poo-2e 2 A a xget(poo)-A ys (poo)ll < {lQE+2A)t A +Ae 2 e~ T 

(27) 

After a short transient time of the order -log-, the ex- 
ponential term can be neglected. Furthermore, note that 
the internal system dynamics £ sys may be time depen- 
dent and thus encode correlations of different components 
of the system in its time dependence. 



IX. 4D TORIC CODE 
A. The 4D Toric code as a stabilizer code 

We will now provide an informal description of the 
4D toric code. For every vertex of an N x N x N x TV 
lattice, there are 6 orientations of faces on which physical 
qubits arc located. Thus, the 6 x TV 4 physical qubits 
are arranged on the 2D faces of a 4D PBC lattice. We 
can now introduce an over-complete set of local stabilizer 
generators for the code, half of which correspond to ID 
edges, the other half corresponding to 3D cubes. For 
each ID edge, there is a tensor product operator Z® 6 , 
the product of Z operators acting on the six 2D faces to 
which this edge belongs. Dual to this, for each 3D cube, 
there is a tensor product operator X® 6 , the product of 
X Pauli operators over the six 2D faces of the cube. Two 
edge and cube stabilizers overlap iff the edge is an edge 
of the cube, and then their overlap will be in exactly two 
faces. Thus all stabilizer generators are seen to commute. 



This expression may be massaged into a form which may 
be more readily bounded. The steps involved include, 
expanding poi according to eq. ( fl7|), t hen expanding 
appearances of poi according to eq. (|22[) and finally in- 
tegrating numerical factors and grouping terms. After 



B. Logical degrees of freedom 

Counting of the remaining degrees of freedom ad- 
ditional to the stabilizer syndrome obtained is not as 
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straightforward as for the 2D toric code, where every syn- 
drome with an even number of anyons was possible. In 
the 4D case, the required condition is that the set of un- 
satisfied stabilizers is only allowed to be a combination 
of closed loops (in the lattice and dual lattice respec- 
tively). However, one can explicitly construct six pairs 
of anticommuting logical operators which commute with 
all stabilizer terms, one pair for each of the six possible 
plane orientations. From each pair, one operator is a full 
plane of X rotations along a full plane wrapping around 
the grid in one of the six possible orientations. The sec- 
ond operator from each pair consists of a dual plane of 
Z operators arranged along the perpendicular plane ori- 
entation. Although analogous to the logical operators on 
the 2D toric code, this image probably stretches our 2D 
or at most 3D imagination. Thus, to obtain an intuition 
about this construction it is convenient to provide formal 
expressions which one may operate with. 



C. 4D PBC lattice notation 

Each vertex of the 4D periodic lattice can be identified 
by a four component vector v = vq,vi,V2,vs € 7L A N . For 
each vertex v, there are four edges e, six faces p and 
four cubes c having the vertex as a lower corner. These 
orientations may be described by four component binary 
vectors 

e,p,c e {(v ,v 1 ,v 2 ,v 3 ) | w,e{0,l}}, (28) 

with edge e, face p, or cube c orientations satisfying the 
additional condition X)t=o w * equal to 1, 2 or 3 respec- 
tively. Each physical qubit can be identified with a tuple 
v, p, where p identifies the plane orientation and v its 
lower side corner. The Z type edge stabilizers E$,e. are 
given by 



(29) 



eCp 



with six participating physical qubits. Finally, the X 
type cube stabilizer are given by 



X. 



V,p ® -^Q?-|-c— p,; 



(30) 



also with six participating physical qubits. 

We will now describe a set of logical operators com- 
muting with all stabilizers which will be used to encode 
information in absence of errors. There is one pair of 
such anticommuting logical operators for each plane ori- 
entation p and they are given by 



N 



X 



nei+me2,p 



Z, 



n ,77i—l 



(31) 

with ei + £2 = p and e% + 84 = p . It is easy to see 
that according to this definition, the two logical opera- 
tor Xp and Zp anticommute, as they coincide only at 



qubit (0,p). One can further verify that such operators 
commute with the complete set of stabilizers. Finally, it 
is not hard to see, that if one assumes the state to be 
in the code subspace (i.e. +1 eigenstate to all stabiliz- 
ers), then any homologically equivalent surfaces results 
in equivalent definition for the operators. 



D. 4D Quantum Toom's rule 

We now define a local update rule which will later be 
used in two ways, first as a dissipation mechanism capa- 
ble of keeping errors from accumulating too badly, second 
as the basic component of an information recovery proce- 
dure permitting removal of all errors to allow information 
read-out. The update rule is analogous to Toom's rule 
for classical information stored in a 2D lattice. While the 
prescription of Toom's rule is to flip a bit if it is different 
to both its two lower side neighbors, the prescription in 
4D will be to X "flip" a qubit if both its neighboring 
lower side Z edge stabilizers are not satisfied, but also to 
Z "flip" a qubit if both its lower side X cube stabilizers 
are not satisfied. This is, a local rotation may be per- 
formed depending on neighboring stabilizer state. This 
is in complete analogy to an interpretation of Toom's rule 
in terms of local stabilizers. One property that permits 
analytic and numerical analysis of such a scheme is the 
decoupling of recovery for X and Z logical operators. 

For each qubit (v,p), we can write the super-operator 
describing the quantum jump implementing the update 
rule as 



^v.pip) ^v,pPv,pPPv,p^v,p + Py.p pPy,p 



X 



(32) 



where P^p is the projector onto the subspace where a 
Z flip should be performed on qubit (v,p) and P^p the 
orthogonal subspace. Assuming p = e\ +82 the projector 
may be defined as 



1 



(33) 



Analogously one may define an update rule TZ§p which 
in a similar way, introduces an X "flip" depending on the 
corresponding projectors P?p in terms of Z type stabi- 
lizers. 



E. Full recovery and error corrected operators 

The superoperators 71-?- and IZi, always commute. 
Only recovery operators of the same kind may lack com- 
mutation when considering neighboring plaquetes. In 
particular, to define a full recovery operation 1Z in terms 
of these local recovery update rules, it is necessary to un- 
ambiguously specify an order of application. Indeed, in 
our simulation code, a sweep through the lattice is taken 
as this order and we observe a good performance in recov- 
ering the originally encoded observables (Fig. [4J. Once 
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is possible to define robust logical observables X EC and 
Z? c such that 



the recovery operation 1Z is unambiguously specified, it corresponding to the Toom like quantum jump super- 

IX Z\ II 

operators 1Z\~ introduced in Eq. 32 We perform 
numerical experiments to determine the relaxation time 
for logical observables (i.e. tv[X EC p(t)} = tr{X L 1Zp(t)}.) 
Evolutions are taken to start in a code state with an un- 
ambiguous X L or Z L logical value and consistency of the 
error corrected logical observables are checked regularly 
in time. The mean time to the first change in the value 
observed for X EC or Z is taken as an estimator of the 
relaxation time. 



\x{Zf c p) = tv(Z^TZp) 
Or more compactly 



tx(X E c p) = tr(X^Up). 



z? c = H (zi) 



and X? c 



K(X±) 



(34) 



(35) 



Thus, error corrected logical observables (super-index 
EC), provide a robust result when evaluated on a state 
with sufficiently few errors and coincide with logical op- 
erators on the error-free subspace. 




0.04 0.06 
Depolarization fraction 



0.1 



FIG. 4. (color online) Recovery probability of an encoded 
observable in the 4D toric code is plotted as a function of 
depolarization probability per qubit. Odd lattices sizes from 
1 to 11 are represented in the different curves and suggest a 
critical depolarization probability of approximately 7.5%. 



G. Numerical considerations 

Evolution under this master equation can be numer- 
ically simulated efficiently for a commuting set of ob- 
servables such as the edge stabilizers and a commuting 
set of logical observables. This means that a classical 
Monte Carlo simulation is enough to study the proba- 
bility of obtaining the correct outcome when measuring 
distinct logical observables which were initially well de- 
fined. The results of such simulations are presented in 
(Fig. [2| for different sizes of the lattice grid up to 11 
and different values for r e . These suggest a critical value 
for the noise rate T* « 0.004, below which arbitrarily 
long relaxation times may be achieved by increasing the 
lattice size. Given that below threshold error rates, the 
information lifetime seems to grow exponentially with the 
lattice size, and that the simulation time per unit time is 
also proportional to the fourth power of the lattice size, 
it is numerically costly to extend our evidence to larger 
lattices. 



H. Definition of efficient recovery 7Z 



F. Master equation 

We study a master equation including a locally depo- 
larizing noise term of strength r c , and the proposed Lind- 
blad terms intended to avoid error clusters from growing. 
The simulated master equation may be written as 

P = £P = FCiDToomP + r e £de P P (36) 
where the dissipative protection CiDToom is given by 

C-iDToomP = /^2 L i!,pP L Z.l - 2 { L Z,i L 0,P> p\ 

v,p - (37) 

+ L lpp L Z-\{ L Z L lpp} + - 

The protecting Lindblad operators are 



To check whether the encoded observable is still recov- 
erable at time t, we apply a correction super-operator 7Z 
on p{t) . The definition of 1Z consists of sequentially ap- 

plying the local jump superoperators 1Z\ p in a sweep- 
ing order. This shows a high performance for removing 
all error domains (i.e. it presents a numerical threshold 
to depolarizing noise on up to ~ 7.5% of the qubits as 
shown in (Fig. [4| .) Furthermore, computer simulations 
of 1Z are efficient, requiring a minimal amount of 0(N 4 ) 
operations. This is important since the 7Zp(t) must be 
checked every unit time to obtain an estimate of the re- 
laxation times of error corrected logical observables. 



X. CONCATENATED-CODE DISSIPATION 

Paz and Zurek [T5] presented the first studies of pro- 
tecting quantum information through the use of a con- 
tinuous dissipative process. They introduce a general 
master equation form for the class of stabilizer QECC 
and analyze their performance in some simple cases. In 
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this section, we will adapt their construction and pro- 
pose master equations for concatenated QECC which are 
provably robust against sufficiently weak local noise. 

For stabilizer QECC, the recovery super-operator 1Z 
can be written in Kraus form as 



H{p) = E R {m) P {m) P P im ^ R {m) \ (39) 

(m) 



where p( m ) are projectors onto orthogonal syndrome 
subspaces with J2( m ) = ^ anc ^ are unitary 

recovery operators of tensor product Pauli form. The 
operators Z/ m ) = p( m )p( m ) can be interpreted as Lind- 
blad operators to give way to a protecting master equa- 
tion. However, this approach can not provide more than 
a constant improvement in the relaxation time for logical 
observables. It can be seen that given an error rate and a 
correction rate, there is an upper bound on the relaxation 
time of the logical encoded bit which is independent of 
the code and the number of physical subsystems it uses. 

We propose extending this master equation model to 
one which allows performing many such recovery opera- 
tions in parallel. In the case of concatenated QECC these 
will correspond to error correction at the different levels 
of concatenation. Recovery operations at the same level 
of concatenation act independently of each other as they 
involve disjoint subsets of physical qubits. Most of the 
work goes into designing recovery operations at different 
levels of concatenation which do not interfere undesirably 
(i.e. commute) and proving that they provide a similar 
protection from local errors to the one achieved by con- 
catenated QECC in the circuit model. 

We will define a dissipative concatenated quantum 
memory based on a [[k, 3, 1]] QECC. An M-level encod- 
ing will thus make use of k M physical qubits. A labeling 
for each physical qubit may be given in terms of an M 
component vector v G Zjj^ (i.e. with each component 
going from 1 to k. Partial vectors v with M — I compo- 
nents will identify mutually disjoint blocks of k l physical 
qubits. Thus if v denotes a particular set of k l physical 
qubit, then the vector vo : v, with one additional compo- 
nent vq and identifies a sub-block of k 1 ^ 1 physical qubits. 
The number of components of a vector v will be denoted 
by | v | , with being the unique zero component vector. 

A stabilizer QECC on k qubits can be characterized by 
the definition of the Stabilizers the projectors onto 
syndrome subspaces P^ and the corresponding error re- 
covery operators the logical operators X L ,Y L , Z L , 
the recovery super-operator 1Z and the error corrected 
Pauli observables X EC , Y EC , Z EC . It is instructive to 
present the definition of these objects for a simple QECC 
making it easier to latter provide the recursive definitions 
required for the concatenated QECCs. These definitions 



are given by 





= sl w <g> s2 (f) <g> . . . <g> sk ij) 


pU) 




RW 


= rl (j) <g> r2 W ® . . . <g> rk {j) 


K{p) 


= Y,R {k) P {k) pP {k) R {k) 


a L 


= crl <g> <t2 ® . . . ® ak 




= TZ(a L ). 



(40) 



The coefficients relating stabilizer operators with 

specific projectors. Lowercase Latin letters as well as cr, 
stand for one of the four single qubit Pauli operators 
{1,X,Y,Z}. Thus er L , is a logical operator on the code 
and as a stabilizer code can be expressed as a tensor 
product of single qubit operators. Finally A denotes the 
super-operator dual to A (i.e. if A(p) = J2k AkpA\, then 

A(0) = £ fe 40A fe ). 

We may now give the analogous definitions for the case 
of an M level concatenated code. Here, objects must be 
further specified by a vector v of at most M components 
indicating the physical qubit or group of qubits they act 
on. Some of these objects require a base case definition 
for | v\= M, 









= O-u 


.EC 


= O-u 


V 


j v 


= 1. 



(41) 



In this case, v identifies on which physical qubit(s) the op- 
erators act on. The super-operator T represents the full 
recovery operation which is trivial in the case of physical 
qubits. For the rest of the objects, definitions are only 
required for | v \< M. 





= (sl (j) )f£® 


...®{sk^) EC , 


P U) 

V 








= (rl«))f :ff ®. 


..®{rk^)lj 




= E4 fe) 4 fe) 


pPi k) Ri k) 

r v V 








or 








= TZ^o {T 1:i j(g 


i ... (8i J"^) . 



The main distinction from non-concatenated definitions 
is that the subindex v has been incorporated everywhere. 
In addition, tensor product decomposition of operators 
now runs either in terms of logical operators (super-index 
L) or error corrected observables (super- index EC). Fi- 
nally, a distinction is made between J 7 ^, which corrects 
all errors in a given block of qubits denoted by v and 
1Z$ which corrects for only the highest level errors within 
that block. This distinction may seem artificial since 
for a simple code (i.e. | v \= M — 1) a full correction 
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corresponds to correcting the highest level error blocks 
possible. 

We will now concentrate on some of the properties 
these recursive definitions carry that will later allow us 
to define the dissipative concatenated QECC and prove 
robustness results. The main property relating logical 
and error corrected operators verified by definition is 



ti[<4°p] = tr[<#F*(p)]. 



(43) 



The meaningfulncss of error corrected operators thus 
stems from the fact that if relatively few errors are ap- 
plied to an encoded state, the error corrected operator 
provides the same expectation value as the logical oper- 
ator on the uncrrcd state 



Error terms £> n oise.ir are single qubit superoperators with 



norm bounded by ||X> no j aei i?|| < IV 
dissipation T> correct ^ is defined by 



The protective 



T^correct,v(p) — r correct, v[R-v{p) P] (50) 

which can be written in Lindblad form as 

v correct A P ) = £4 j) /4 j)t \ {4 j)t 4M + (si) 

0) 



with Lindblad operators 

V correct,v I *'tf r 



(52) 



tv[0% c £(p)} = tr[0±p] Mp e codespace^, 



(44) 



provided that the error super-operator £ contains only 
"few error" Kraus operators. More precisely, the expec- 
tation values are equal provided that the Kraus opera- 
tors for £ contain less than ^^ilJ-M-M errors. More can 
be said in terms of the structure of correctable errors. 
Namely, there is a constant error threshold provided a 
random distribution of uncorrelated errors is assumed. 

Another key property which can be guaranteed induc- 
tively is that the commutation/anticommutation relation 
between logical operators and error corrected observables 
should be the same as between bare operators. 



[<rlf c ,<72£]±=0 [al,a2] z 







(45) 



An even stronger statement can be made about products 
of logical operators (error corrected observables) 



erlcr2 = OaS 



alia2i = QaZh and 
CT lf <V 2 f c = 6a3§ c 



(46) 



where 9 is a phase in {1,-1, i, — i}. 

The projector operators at each level are related to the 
presence of logical errors at the immediately preceding 
level. This can be seen through the identity 



i?S 0) = 1. 



(47) 



which relates P^\ the trivial syndrome projector to 
other syndrome projectors. The relation of this projector 
with the recovery operations is captured by 



P^n s (p)p^ 



(48) 



The master equation. - considered for a dissipative pro- 
tection on a concatenated QECC will contain error terms 
T> no ise,v on single physical qubits v as well as correction 
terms corresponding to each of the blocks. The full mas- 
ter equation reads 



p= Yl v - 

\v\=M 



MP) + V correct ,v{p)- ( 49 ) 



|v|<M 



We will prove the robustness of the highest level 
observables X^ c ,Y^ C 7 Z^ c under the combination of 
weak local noise and this dissipative protection. To do 
this, we focus on the observables {Pjf* ■ \v\ < M}. To- 
gether with an arbitrary error corrected observable at the 
highest level, these constitute a complete set of quantum 
numbers. The most attractive features of these observ- 
ables is that both single qubit Pauli errors and the re- 
covery operations may be described by classical deter- 
ministic transition rules in terms of this specific set of 
quantum numbers. Furthermore, the events influencing 
each of these quantum numbers may be simply charac- 
terized. Namely, only recovery or physical error events 
located at«^? can influence the validity of ■ This 
will allow us to provide upper bounds for the probability 
of introducing logical errors. 

It is useful to define certain additional projectors in 
terms of the set of commuting projectors {Pi^ : M > 
1*1} 

HasError(tf) = 1 - pi 0) 

IsError(j : v) = pi X]) + P { ? ]) + pi Zj) (53) 
Enablcd(j : v) = 1 - pi a) - pi X]) - pi Yj) - pi Zj) 



The recovery operation 1Z$ has a non trivial effect only 
for the subspace "HasError(w)" . Furthermore, in the sub- 
space "IsError(j : v)" , the effect of applying the recovery 
operation 1Z$ is to apply a logical operation on j : v. 
Finally, the projector "Enabled(j : v)" the difference be- 
tween the two and indicates that there is already a logical 
error among the immediate components of v, but that it 
is not at j : v. This last projector will be instrumental in 
bounding the probability for physical errors to be raised 
as logical errors. In the case of the perfect five qubit 
code, it is a necessary and sufficient condition for a logi- 
cal operation at j : v be seen (in terms of the stabilizers) 
as raising a logical operation at v. The following short 
hand notation will be used to express the probability of 
satisfying these predicates (projectors) 



(P) t =tr[Pp(t)]. 



(54) 
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A. Bounding error probabilities 

of constitute the core of proving the robustness of error 
corrected observables under such a dissipative dynamics 
as Eq. (49 1. In particular, we wish to prove inductively 
that 

Vt (HasError (u)^ < p n where n = M — \v\. (55) 

Since the initial state p(0) is by Hypothesis a code state 
at all levels, we have that 



VtT : (HasError(i7)) f= o = 0. 



(56) 



The trick now is to obtain an upper bound on the rate at 
which these probabilities may increase and upper-bound 
the actual probability by a fixed-point value. Let us first 
illustrate this method by considering a simple example 
provided by \v\ = M — 1. 



< 



d(HasError(z7)) t 
dt 

r c orrcct,{?(HasError(if))t 



(57) 



Note that we have excluded processes by which a physical 
error cancels a preexisting error. 

From the rate bound, we may extract a fixed-point 
upper-bound and use it to bound the actual probability 



(HasError(i7)) t < 



fcrv 



fcr n 



(58) 



Assuming r correct ,s > r correct)Af _|^ , we may further sim- 
plify the bound to 



(HasError(i?)) t < 



--■■Pi- 



(59) 



correct, 1 

We may take a similar approach to bound the rate at 
which errors accumulate at higher levels (i.e. M — \v\ — 
n + 1). However, the expressions required here are a bit 
more complicated. 

d(HasError(u)) t 
dt 



(60) 
(61) 



< Y r "°isc,«( II Enablcd(w))j 

\u\=M 1 

- Incorrect , v (HasError ( v) ) t 

<r noisc Y II (Enabled(w;)) t (62) 

\u\=M 

- r c orrcct , v (HasError [v) ) t 

<r noisc II (HasError(w;)) t (63) 

\u\=M 

- r corrcct (HasError (v)) t 

n 

< k n+1 T noise Y[ Pj (64) 
-r c0 rrcct,n+i (HasError (v)) 



A non trivial step is taken in going from (Ml to 62 , 



where the probability of a conjunction is taken to be a 
product of probabilities (i.e. independent probabilities) . 
This property will be proven in appendix |XI| 
In turn, this leads to the fixed point bound 

p+ir FT™ n 
(HasError(w)) t < J - =: p n +i- (65) 

1 correct, n+1 

From here, we inductively derive the expression 



Pn 



x noise ^ 



2"-l 



correct, n llj = l 1 correct,;/' 



Making the additional assumption r correc t j:/ - 
we may simplify this expression to obtain 

/ 1 noise ^ ^ 

Pn 



= r c 



,6 2 



5_ 



(66) 



(67) 



In turn, this tells us that if r no j sc < (<5/fc) 2 r corrcc t, then 
the probability of having non trivial syndrome decreases 
double exponentially with the level of the syndrome. 

Our final goal is to obtain an expression bounding the 
rate at which logical errors are introduced. One possibil- 
ity, is to study the decay rate for any of the three highest 
level logical Pauli observables. Since these three consti- 
tute a full set of observables for the logical subsystem, 
their preservation implies high fidelity storage of quan- 
tum information |33j . 

A logical error or flip of the highest level logical observ- 
ables, can be introduced whenever a physical error occurs 
at a site which is enabled to raise the error at all levels. 



Employing bounds similar to those in Eqs. (61 1-(64) one 
arrives at 



dt 



< 



\u\=M 



noise ,ii 



Yl Enabled (w) 



\uy-wyi 



r +f 5 2 

x correct u 



(68) 



(69) 



(70) 



indicating that for a sufficiently low physical error rate, 
the logical error rate is suppressed double exponentially 
in terms of M, similar to results for concatenated QECC 
in a quantum circuit model. 



XI. PROOF OF INDEPENDENCE FOR THE 
ENABLED PROPERTY 



In this section we assume a Pauli noise model and prove 
that the Enabled property along the different truncations 
of the same physical address are statistically indepen- 
dent. More specifically, the factorization 

/ Yl Enabled(w)) = JJ (Enabled^)) t (71) 
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holds for a noise process of Pauli form 

T^noise,v{P) = X! (VvPVv ~ P) ■ (72) 

cr£{X,y,Z} 

The restriction of the noise process to Pauli form Eq. 



( 72 ) is clearly undesired. However, it provides a sufficient 



condition to prove Eq. ( 71 1 , which does not hold for 
general noise. We expect the need for this assumption 
to be an artifact of our proof technique and that our 
main result, i.e. Eq. (68), will essentially hold for any 



independent noise model. 

The proof relies on the independence of the different 
processes which introduce physical errors and perform 
recovery operations. An event Ev^ will be associated 
to each vector w, with \w\ — M corresponding to the 
introduction of physical errors at w and < M cor- 
responding to recovery operation TZ^;. Each event Ev^ 
can be seen as the state dependent application of a ten- 
sor product Pauli operator. Furthermore, for \w\ < M 

(i) 

the operator only depends on the quantum numbers P^j 
and must be a logical Pauli operators at some wq : w. 
The correction operators satisfy this property by design. 
In turn, for \w\ — M, Ev^j applies a randomly chosen 
physical Pauli operator at w according to the Pauli form 
noise model Eq. (72). It can be seen that under these 



conditions, only events Ev^ such that w y v can directly 

(i) 

affect the quantum numbers P$ . Thus, given a history 
L of events Ev^j applied to an initially encoded state, the 
quantum numbers P$ are well defined and depend only 
on the sub-history of events L' containing the events Ev^ 
with w )p v . 

Since Enabled(i>o : v) can be defined in terms of the 

(i) 

py it may only depend on the sub- history of events 
Ev^j with w )p c. Furthermore, Enabled(i>o : v) will be 
shown not to depend direct or indirectly on events Ev^ 



with u )p vq : v. This can be seen as a consequence 
of Enabled(i>o : v) commuting with any Pauli operator 
acting on qubits w with w )p Vq : V. 

Proving Eq. (TTTT) may be split in the following steps 



( | | r J nablea(u7j ) 


(73) 


\ A / j \ j r 1 T T7> ii j / -*\ r 1 

> PL{t) ti[ (1 _Lnabled(u;)Lp J 


(74) 


^Pi(t) Jl tr[Enabled(u?)Lp ] 


(75) 


^2Ph{t) J| tr[Enabled(w)L^p ] 


(76) 


11 5Z PL ™ W tr [ Enab led(wJ)i^p ] 


(77) 


Yl (Enabled(w)) t , 


(78) 



which will be subsequently explained and justified. As 
a first step, the master equation defining p(t) is unrav- 
eled [24] into event histories L to obtain Exp. (74). 



Given that every event history L implements a Pauli op- 
erator which produces eigenstates to all the projectors 
Enabled(w), the 0,1 expectation values may be factor- 
ized to obtain Exp. (75). Expectation values depend 
only on disjoint sub-histories L Wo .$ ( a history of events 
uniquely determined by filtering events Ev^ such that 
u y w but not u >~ wq : w from L ), leading to Exp. 
(76). Furthermore, the sum of p^ consistent with given 



sub-histories may be written as a product of the in- 
dependent probabilities pl^ of such sub-histories, thus 
leading to Exp. (77). Finally, each factors in Exp. (76) 



may be seen to be the history unraveling of each of the 
factors in Exp. (78), which is what we set out to prove. 



